Diffractive 0-meson photoproduction on the proton near threshold 
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Photoproduction of 0-meson on protons was studied by means of linearly polarized photons at 
forward angles in the low-energy region from threshold to E-,= 2.37 GeV. The differential cross 
sections at t = — |t| m i n do not increase smoothly as _E 7 increases, but show a local maximum at 
around 2.0 GeV. The angular distributions demonstrate that </>-mesons are photo-produced predom- 
inantly by helicity-conserving processes, and the local maximum is not likely due to unnatural-parity 
processes. 

PACS numbers: 13.60.Le, 25.20.Lj 



The gluonic aspect of quantum chromodynamics 
(QCD), especially glueballs, has been of wide interest 
in hadron physics. It has been suggested that there 
is a connection between the glueball Regge trajectory 
(J PC — 2 ++ ,4 ++ ...) and the Pomeron trajectory [1|. 
The diffractive photoproduction of 0-mesons has tradi- 
tionally been used to study the Pomeron exchange pro- 
cess |2| . This is because the baryon and meson exchange 
amplitudes in the s- and t-channels are suppressed by 
the Okubo-Zweig-Iizuka (OZI) rule. Photoproduction of 
</>-meson is useful to study not only Pomeron exchange, 
but also other hadronic interactions mediated by multi- 
gluon exchanges which are difficult to identify in the 
other hadronic reactions due to large contributions from 



baryon and meson exchanges. 

The low-energy diffractive photoproduction of <p- 
mesons was suggested 0, 0] to be sensitive to a daughter 
Pomeron trajectory associated with a glueball (J PC — 
++ ) 0, El H El- Its contribution is expected to de- 
crease rapidly with an increase of photon energy, whereas 
the contribution from Pomeron exchange increases. This 
difference may lead to a non-monotonic energy depen- 
dence of the forward-angle cross section near threshold 
(Ey = 1.57 GeV). The existin g cr oss section data in the 
low-energy region 0, 0, llftllll Il2| are still too poor 
to ascertain a possible signature of such a non-monotonic 
behavior. 

The background contributions from s- and u-channels 
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diagrams, such as direct <fr radiation from the nu- 

cleon |4 1 1 311 and production in nucleon resonance de- 
cay [14], [15J, are predicted to be small at small \t\ 
{t = {p<t> — Pj) 2 )- However, the contributions from the t- 
channel exchanges of pseudoscalar mesons (n, y), scalar 
mesons (/o, ao) 0, 0] and a tensor j' 2 meson ^t) are 
predicted not to be negligible. The energy dependence of 
those meson-exchange processes is expected to be similar 
to that of the daughter Pomeron trajectory. Therefore, 
to determine the relative contributions of these processes, 
we need to analyze the spin observables using linearly po- 
larized photons. 

The spin observables are studied via the decay angular 
distribution of the </>-meson in the K + K ~ decay mode. 
The decay angular distribution W (cos 6, </>,$) is a func- 
tion of the spin density matrix elements |l8| . where 9 
and cf> denote the polar and azimuthal angles of the K + 
in the 0-meson rest frame. The azimuthal angle of the 
photon polarization in the center-of-mass frame is de- 
noted by The relative contribution of natural-parity 
exchange and unnatural-parity exchange is related to 
the density matrix element p\_ 1 (= l/2(p\_ 1 — \mp1_-J) 
which is extracted from the one-dimensional distribution 
W((j> - $) 15] through 

W{4> - *) = -^(1 + 2P 7/ 5l_x cos 2(0 - $)), (1) 

where P 7 is the polarization of the photon beam. The 
available data at E 1 = 2.8,4.7,9.3 GeV and 20- 
40 GeV 0| support the dominance of the helicity- 
conserving natural-parity exchange processes. However, 
there is no measurement of polarization observables near 
the threshold. 

The decay angular distributions also provide infor- 
mation on helicity non-conserving processes. Recent 
measurements at low energies with unpolarized pho- 
tons suggest significant contributions from the helicity 
non-conserving processes at large momentum transfer 
t |l2ll20| . The contribution from helicity non-conserving 
mechanisms is examined by a deviation from the sin 2 9 
behavior in the one dimensional distribution W{cos9) 
and an oscillation in the one dimensional distributions 
W((f> + <P),W((f>) and W($). 

In this Letter, we report measurements of the differen- 
tial cross sections (da/dt) at small \t\ and the first mea- 
surements of decay angular distributions near threshold 
with linearly polarized photons. Linearly polarized pho- 
tons were produced by means of the backward-Compton 
scattering of laser photons off the 8 GeV electron at 
the SPring-8 BL33LEP beamline (LEPS: Laser Electron 
Photons at SPring-8 facility) [2lJ . The maximum energy 
of the photon beam was 2.4 GeV. The photon energy was 
determined by measuring recoil electrons using a tagging 
counter with a resolution (a) of 15 MeV. The typical 
photon flux was about 10 6 s , which was monitored 
by counting scattered electrons with the tagging system. 



The systematic uncertainty in the photon flux measure- 
ment was estimated to be 3%. The degree of linear po- 
larization varied with photon energy; it was 95 % at the 
maximum energy, 60 % at 1.57 GeV. A liquid hydrogen 
target with a length of 50 mm was used in the experi- 
ment. A similar experiment with nuclear targets and the 
associated analysis have been reported in Ref. |2^ . 

The momenta and the time-of-flight (TOF) of pro- 
duced charged particles were measured with a magnetic 
spectrometer [21| . The angular coverage of the spectrom- 
eter is about ±0.4 rad and ±0.2 rad in the horizontal and 
vertical directions, respectively. The momentum resolu- 
tion (a) for 1 GeV/c particles was 6 McV/c. The TOF 
resolution (tr) was 150 psec for a typical flight path length 
of 4 m. The mass resolution (a) was 30 MeV/c 2 for a 
1 GeV/c 2 kaon. Pions with momenta higher than 0.6 
GeV/c and e + e~ pairs were rejected by using an aerogel 
Cherenkov counter in the trigger level. An over-veto rate 
in the trigger was estimated to be less than 2.1%. 

The incident photon energy and the momenta of 
K + K~ tracks or K p tracks were measured to identify 
the reaction — > K + K~p followed by the <j> — > K + K~ 
decay. Based on the detected particles, we define two 
types of event topology: A + A"~-reconstructed events 
(KK mode) , and A^p- reconstructed events ( Kp mode) . 

The missing mass distribution for the p("f, K + K~)X 
reaction (denoted as MM(j,K + K~)) is shown in 
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FIG. 1: (a) Missing mass distribution for the p("f, K + K~)X 
reaction in KK mode, (b) Missing mass distribution for the 
p(7, K ± p)X reaction in Kp mode, (c) and (d) are the K + K~ 
invariant mass distributions after the cut on the missing mass 
for KK and Kp modes, respectively. The hatched histograms 
are the simulated background. 
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FIG. 2: Differential cross sections for the jp — ► (j>P reaction. The dashed curves are the results of the fit using an exponential 



function {{da/dt) 



t=—\t\. 



Ht+\t\ min ) 



) with {d<7 /dt)t = -\t\ m n and b as free parameters. The solid curves are fitted results with 



fixing b = 3.38 GeV . The error bars represent statistical errors. The systematic errors are discussed in the text. 



Fig. nja) for KK mode. A sharp peak corresponding 
to the proton was observed with an average mass resolu- 
tion (<r) of 10 MeV/c 2 . The missing mass distribution 
for the p(j, K ± p)X reaction {MM{j, K^p)) is shown 
in Fig. ^b) for Kp mode. A clear peak at the kaon 
mass was observed with an average mass resolution {a) 
of 10 MeV/c 2 . In Kp mode, there were contributions 
from non-K + K~p final states. The background below 
the kaon peak is attributed mainly to hyperon photopro- 
duction having a non-K + K~p final state, such as Kpn^y 
, Kpinr. The background above the kaon peak is due to 
KKpir events. A 3 a cut on the missing mass spectrum 
was applied to select the K + K~p final state. 

Figure Qfc) and (d) show the K + K~ invariant mass 
distributions for KK and Kp modes, respectively. In Kp 
mode, the momentum of the missing kaon was calculated 
by assuming a K + K~p final state. The cut point on the 
K + K~ invariant mass was set to 1.009 < M(K + K~) < 
1.029 GeV/c 2 which corresponded to about 10 % loss of 
<f> events. The background in the <f> peak region was es- 
timated with the following method. We considered two 
sources of background; photoproduction of A(1520) and 
a K + K~p final state without forming any narrow reso- 
nance structure in either K + K~ or K^p system (non- 
resonant KKp). The background level was estimated 
from the yields below and above the c6-meson peak by 
using Monte Carlo simulations which were fitted to the 
angular distributions of K + , K~ and p in the real data. 
The Monte Carlo simulations reproduced the K + K~ in- 
variant mass (Fig.^c) and (d)) and K~p invariant mass 
distributions in the real data. Although there is a kinc- 
matical overlap of A(1520) and <j> production in this en- 
ergy range, the contamination of A(1520) events in the 
final sample is suppressed in small \t\ regions. It was 
estimated as less than a few percent. The estimated sys- 
tematic error on the cross section due to the background 
subtraction procedure was less than 0.8%. 



The acceptance of the spectrometer was determined in 
Monte Carlo simulations using the GEANT3 simulation 
package j^l ■ Geometrical acceptance, resolution and ef- 
ficiency of the detectors were taken into account. Since 
the acceptance depends on the input distributions, the 
simulations were iterated to reflect measured da/dt and 
angular distributions, having started from flat distribu- 
tions. The acceptance also depends on beam polarization 
and the mode of reconstruction. The validity of the ac- 
ceptance calculation and the background subtraction was 
confirmed by checking the consistency of the cross sec- 
tion results among different reconstruction modes, and 
also by checking the consistency of the decay angular 
distributions obtained with different beam polarizations. 

The differential cross sections were measured in terms 
of t + |f | m j n where |f | m j n is the minimum 4-momentum 
transfer from the incident photon to the ^>-meson. Fig- 
ure[21shows the da/dt in different photon energy regions. 
The da/dt showed a forward peaking shape, suggesting 
the dominance of t-channel exchange processes. A fit to 
the da/dt was performed with an exponential function; 
i.e. {da/dt) t= „ Wmm e b( - t+ ^™"^ with {da / dt) t= _\ t \ mm and 
b as free parameters. No strong energy dependence of the 
slope b was found beyond statistical errors. The average 
value of the slope b was 3.38 ± 0.23 GeV" 2 . When the 
average slope for data at all energies was used in the fit, 
the fitting curves described the data points well. 

Figure [21 shows the energy dependence of 
(da/dt) t —_i t i in when b is set to the average slope. 
The energy dependence of {da / dt) t= \ t \. n shows a 
non-monotonic behavior with a local maximum at 
E 1 ~ 2 GeV. The local maximum is also seen using only 
the differential cross sections at the lowest t + \t\ m i n 
bin where acceptance and signal-to-noise ratio is at a 
maximum. The local maximum still persisted when the 
analysis was repeated excluding events near the A(1520) 
peak in K~p system. 
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FIG. 3: Energy dependence of (dcr/dt) t= _i t \. . The closed 
circles are the results of the present work. Other data points 
are taken from Ref. H H 113 [H [13 ■ The error bars rep- 
resent statistical errors. The systematic errors are discussed 
in the text. The solid curve represents the prediction of a 
model including the Pomeron trajectory, n and r\ exchange 
processes fl5||. 



local maximum of the cross section 
y < 2.17 GeV) and (2) above the local 



Data were compared with the prediction of a model 
including Pomeron exchange, ir and rj exchange pro- 
cesses [lq. A x 2 test was performed to check whether 
the model prediction was statistically compatible. It gave 
X 2 = 140 for 8 degrees of freedom using the present mea- 
surements. The model is inconsistent with the present 
data points although it describes the data rather well at 
higher energies. 

The decay angular distributions in the Gottfried- 
Jackson frame were obtained at forward angles (—0.2 < 
t + \t\min < GeV 2 ) in two different energy-regions: 
(1) around the 
(A£i:1.97 < R 
maximum (Ai?2:2.17 < E 1 < 2.37 GeV) where there are 
enough statistics and the acceptance is fairly flat over all 
angular variables. 

Figure ^a) shows the angular distribution W{aosd). 
In both energy regions, W(cos0) behaves as ~ 
(3/4) sin 2 0, indicating the dominance of helicity- 
conserving processes. A contribution from tensor-meson 
exchange, such as j' 2 -meson exchange, must be small 
since a contribution of this term would result in a de- 
viation from the sin 2 6 form 0| . Note that this result 
is different from the measurement in a wider t range [l^ 
which shows strong violation from sin 2 9 form. This may 
be understood by the production mechanism discussed 
in Ref. El. 
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FIG. 4: Decay angular distributions for -0.2< t+ \t\ m in in the 
Gottfried- Jackson frame. The solid curves are the fit to the 
data. The hatched histograms are systematic errors. 



Figure^fb) shows the distribution W{4>— $). We found 
0.197 ± 0.030(stat.) ± 0.022(sys.) in A£?i and 
0.189 ± 0.024(stat.) ± 0.006(sys.) in AE 2 . The positive 
value for p\_ 1 indicates that the contributions from nat- 
ural parity exchange are bigger than those for unnatural 
parity exchange (7r,?y-meson exchange). The p\_i is the 
same in the two energy regions within errors. This im- 
plies that the relative contribution of natural parity ex- 
change and unnatural parity exchange remains constant 
in the two energy regions. Therefore, it is difficult to 
attribute the origin of the local maximum in the cross 
section to different strengths of the unnatural-parity ex- 
change processes in the two energy regions. 

Other one dimensional angular distributions W (</)), 
W{<j) + $) and W($) are depicted in Fig. He). No 
strong oscillation was found, except that the distribu- 
tion W(<f) at AEi bin showed an oscillation (p1_ 1 — 
0.120±0.027(stat.)±0.011(sys.)). p° 1 _ 1 reflects the double 
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spin-flip transition from the incident photon to the out- 
going 0-meson [l^. The spin-flip amplitudes are exactly 
zero in the case of pure scalar meson exchange and pseu- 
doscalar meson exchange processes. The oscillation in the 
W(<fi) distribution might be understood in the framework 
of a modified Donnachie-Landshoff Pomcron model mo- 
tivated by the non-perturbative two-gluon-exchange dy- 
namics [13 • However, this model fails to reproduce the 
non-monotonic energy dependence (see the solid curve in 
Fig.EJ. 

An alternative explanation might be the manifestation 
of a daughter Pomeron trajectory. In this case, the de- 
cay angular distributions may be similar to those for the 
Pomeron trajectory as observed, since contributions from 
both trajectories involve exchanges of natural-parity par- 
ticles. While the decay angular distributions are useful 
to discriminate natural parity exchange from unnatural 
parity exchange, they are not useful to disentangle the 
two possible natural parity exchanges, i.e. Pomeron ex- 
change and a daughter Pomeron exchange. On the other 
hand, the energy dependence of the cross sections is a 
good indicator for a daughter Pomeron exchange process. 
However, the fit suggested in Ref. Q failed to predict the 
local maximum in the cross section with the proposed set 
of parameters. 

In summary, the photoproduction of the </>-meson was 
studied for the first time by means of linearly polar- 
ized photons at forward angles in the low energy region 
from the threshold energy of # 7 =1.57 GeV to 2.37 GeV. 
The differential cross sections at t = — \t\ m i n go non- 
mo notonically as a function of E 1 , and show a local max- 
imum at around 2.0 GeV. The polar angle distributions 
demonstrate dominance of helicity conserving processes 
and disfavor tensor meson exchange. The azimuthal 
angle distributions over the local maximum suggest that 
the local maximum is not due to additional unnatural- 
parity processes, but likely due to new dynamics which 
may involve a multi-gluon exchange beyond the Pomeron 
exchange process. Further theoretical and experimental 
studies are of great interest for clarifying photoproduc- 
tion mechanisms in the low-energy region with the local 
maximum. 

The authors thank the SPring-8 staff for support- 
ing the BL33LEP beam line and the LEPS experiment. 
We thank H. Toki and A. Hosaka (RCNP) for fruit- 
ful discussions. This research was supported in part by 
the Ministry of Education, Science, Sports and Culture 
of Japan, by the National Science Council of Republic 
of China (Taiwan), Korea Research Foundation(KRF) 
Grant(2003-015-C00130) and National Science Founda- 
tion (NSF Award PHY-0244999). 



* Present address: Department of Physics and Astronomy, 
Ohio University, Athens, Ohio 45701 

' Present address: Department of Physics, Kyoto Univer- 
sity, Kyoto 606-8502, Japan 

* Present address: Laboratory of Nuclear Science, Tohoku 
University, Sendai, Miyagi 982-0826, Japan 

§ Present address: Center for Nuclear Study, University of 

Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 
1 Present address: Department of Applied Physics, Na- 
tional Defense Academy, Yokosuka 239-8686, Japan 

** Present address: Department of Physics, Tokyo Institute 
of Technology, Tokyo 152-8551, Japan 

" Present address: Department of Physics and Astron- 
omy, Rice University, 6100 Main St. Houston MS 108, 
TX 77005-1892, USA 

** Present address: Department of Physics, Tohoku Univer- 
sity, Sendai, Miyagi 980-8578, Japan 

^ Present address: Joint Institute for Nuclear Research, 
141980, Dubna, Russia 

" ■ Present address: Akita Research Institute of Brain and 

Blood Vessels, Akita 010-0874, Japan 
*** Present address: Institute of Statistical Science, 

Academia Sinica, Nankang, 115 Taipei, Taiwan 
' ' ' Present address: Illinois Institute of Technology, 

Chicago, Illinois 60616, USA 
* * * Present address: Research Center for Nuclear Physics, 
Osaka University, Ibaraki, Osaka 567-0047, Japan 

§§§ p resen t address: National Superconducting Cyclotron 
Laboratory, Michigan State University, East lansing, MI 
48824-1321, USA 

[1] S. Donnachie, G. Dosch, P. Landshoff, and O. Nacht- 
mann, Pomeron Physics and QCD (Cambridge Univer- 
sity Press, 2002), reference therein. 

[2] T. H. Bauer, R. D. Spital, D. R. Yennie, and F. M. Pip- 
kin, Rev. Mod. Phys. 50, 261 (1978). 

[3] T. Nakano and H. Toki, in proceedings of International 
Workshop on Exciting Physics with New Accelerators Fa- 
cilities, SPring-8, Hyogo (World Scientific) p. 48 (1997). 

[4] A. I. Titov, T. S. H. Lee, H. Toki, and O. Streltsova, 
Phys. Rev. C 60, 035205 (1999). 

[5] L. S. Kisslinger and W.-H. Ma, Phys. Lett. B 485, 367 

(2000) . 

[6] F. J. Llanes-Estrada, S. R. Cotanch, P. J. de A. Bicudo, 
J. E. F. T. Ribeiro, and A. P. Szczepaniak, Nucl. Phys. 
A 710, 45 (2002). 

[7] J. Ballam et ai, Phys. Rev. D 7, 3150 (1973). 

[8] H. J. Besch et al, Nucl. Phys. B 70, 257 (1974). 

[9] H. J. Behrend et ai, Nucl. Phys. B 144, 22 (1978). 
[10] D. P. Barber et ai, Zeit. Phys. C 12, 1 (1982). 
[11] E. Anciant et ai, Phys. Rev. Lett. 85, 4682 (2000). 
[12] J. Barth et ai, Eur. Phys. J. A 17, 269 (2003). 
[13] Y.-S. Oh and H. C. Bhang, Phys. Rev. C 64, 055207 

(2001) . 

[14] Q. Zhao, B. Saghai, and J. S. Al-Khalili, Phys. Lett. B 

509, 231 (2001). 
[15] A. I. Titov and T.-S. H. Lee, Phys. Rev. C 67, 065205 

(2003). 

[16] R. A. Williams, Phys. Rev. C 57, 223 (1998). 
[17] J. M. Laget, Phys. Lett. B 489, 313 (2000). 
[18] K. Schilling, P. Seyboth, and G. E. Wolf, Nucl. Phys. B 
15, 397 (1970). 



6 



[19] M. Atkinson et al., Z. Phys. C27, 233 (1985). 
[20] K. McCormick et ai, Phys. Rev. C 69, 032203 (2004). 
[21] T. Nakano et al, Nucl. Phys. A 684, 71 (2001). 
[22] T. Ishikawa et al, Phys. Lett. B 608, 215 (2005). 
[23] R. Brun et al, CERN Program Library Long Writeup 
W5013, CERN Applications Software Group (1993). 



